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ABSTRACT 

The existence of a primordial magnetic field at the last scattering surface may 
induce a measurable Faraday rotation in the polarization of the cosmic microwave 
background. We calculate the magnitude of this effect by evolving the radiative 
transfer equations for the microwave background polarization through the epoch 
of last scatter, in the presence of a magnetic field. For a primordial field amplitude 
corresponding to a present value of 10 _9 G (which would account for the observed 
galactic field if it were frozen in the pre-galactic plasma), we find a rotation angle 
of around 1° at a frequency of 30 GHz. The statistical detection of this signal is 
feasible with future maps of the microwave background. 
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1. Introduction 

The origin of the observed /zG magnetic fields in spiral galaxies has been a long-standing 
puzzle for more than three decades (Rees 1987; Kronberg 1994). On the one hand, these 
magnetic fields could have resulted from an exponential dynamo amplification of a small 
seed field with an e-fold period of a galactic rotation (Parker 1979; Zel'dovich, Ruzmaikin, 
& Sokoloff 1983; Field 1994). Alternatively, these fields may have already been frozen in the 
primordial plasma before galaxies formed (Hoyle 1958; Piddington 1964, 1972; Ohki et al. 
1964), and could have even affected the evolution of structure in the universe (Wasserman 
1978; Kim, Olinto, & Rosner 1994). The origin of the primordial field is still a subject of 
speculation (see, e.g. Turner & Widrow 1988; Quashnock, Loeb, & Spergel 1989; Vachaspati 
1991; Ratra 1992). 

In the past, various indirect theoretical arguments were used to favor the dynamo 
amplification mechanism over the primordial origin alternative (Zel'dovich, Ruzmaikin, & 
Sokoloff 1983). However, recent theoretical studies argue that a galactic dynamo should 
saturate due to the rapid growth of a fluctuating small-scale field before it can actually result 
in a coherent large-scale field of the type observed in galactic disks (Kulsrud & Anderson 
1992; Vainshtein & Cattaneo 1992; Vainshtein, Parker, & Rosner 1993; Cattaneo 1994). 
The view that the galactic field may, in fact, be primordial gains additional support from 
observations of damped Lja absorption systems in QSO spectra at ^ abs ~ 2. These systems, 
which are thought to be the progenitors of galactic disks (Lanzetta et al. 1995), add Faraday 
rotation to their background QSOs, consistent with them having /iG fields (Welter, Perry, 
& Kronberg 1984; Wolfe, Lanzetta, & Oren 1992; Kronberg, Perry, & Zukowski 1992). 
Since these systems exist at an epoch when their rotation period is not exceedingly small 
compared to their age, this result could pose a problem to the standard dynamo hypothesis. 
However, the current data may not be sufficient to draw any firm statistical conclusions 
about absorption systems beyond a redshift of 0.4 (Perry, Watson, & Kronberg 1993). 

The potential existence of a primordial magnetic field is also consistent with observations 
of clusters of galaxies. Faraday rotation measurements of radio sources inside and behind 
clusters indicate strong magnetic fields in many of them (Kim et al. 1990, 1991; Taylor 
& Perley 1993). The detected cluster fields have a typical magnitude of a few fiG and a 
coherence length of 10 1-2 kpc. The cores of several clusters contain tangled magnetic fields 
with amplitudes as high as ~ 10 1_2 /xG (Dreher et al. 1987; Perley & Taylor 1991; Taylor & 
Perley 1993; Ge & Owen 1993). In the outer halos of clusters, lower limits ^ 0.1/iG were 
set on the field amplitude, by combining measurements of synchrotron radio-emission from 
relativistic electrons in these halos together with lower limits on the associated hard X-ray 
emission due to Comptonization of the microwave background (Rephaeli 1979; Rephaeli et 
al. 1987; Rephaeli, & Gruber 1988). 
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If primordial in origin, the /iG galactic field could have resulted directly from the 
adiabatic compression of a cosmological field, B ~ (10~ 10 -10~ 9 )G. Using a sample of 309 
galaxies and quasars with a small intrinsic rotation measure, Vallee (1990) was able to set 
an upper limit of 10 _9 G x (^igmVO.OI) -1 on the magnitude of a cosmological magnetic field 
which is coherent on the scale of the horizon; here, fliGM is the ratio between the ionized gas 
density in the intergalactic medium and the critical density, and h is the Hubble constant H 
in units of 100 km s _1 Mpc -1 . However, this limit is based on Faraday-rotation measures, 
for which the contributions of field reversals along the line of sight average-out in a tangled 
field configuration. The limit is weakened to a value ^ 3 x 10~ 8 G x (Oigm^/O.OI) -1 if the 
cosmic field is coherent only on scales <^ lOMpc (see also Kronberg 1994). 

In this paper, we propose a direct empirical probe of primordial magnetic fields. It is by 
now established that the microwave background should have acquired a measurable level of 
polarization at decoupling (see, e.g. Kosowsky 1996; Coulson, Crittenden, & Turok 1995). 
If a primordial magnetic field is present at the last scattering surface of the microwave 
photons, it will cause Faraday rotation of the direction of linear polarization. Since the 
rotation angle depends on wavelength, it is possible to infer this effect by comparing the 
polarization vector of the microwave sky in a given direction at two different frequencies. 
Since both the magnetic field amplitude (oc [1 + z] 2 ) and the baryonic density (oc [1 + z] 3 ) 
increase rapidly with redshift, this effect could potentially have measurable consequences at 
the high redshift of decoupling. 

We can roughly estimate the expected rotation angle of the microwave background 
polarization as follows. Monochromatic radiation of frequency v passing through a plasma 
in the presence of a magnetic field B along the propagation direction q will have its linear 
polarization vector rotated at the rate 



where e and m are the electron charge and mass, n e is the total number density of 
electrons, and x e is the ionization fraction. (Throughout the paper we use natural units 
with h = c = G = 1.) The optical depth for scattering is of order unity out to the redshift of 
decoupling when polarization is generated, i.e. we may substitute / x e n e dt « 1/or where ot 
is the Thomson cross-section. The rms rotation angle can then be easily estimated from the 
time-integral of equation (|1|) by noticing that Bjv 2 is time-independent and by averaging 
ip 2 oc (B • q) 2 over all possible orientations of B, 







where Bq is the current amplitude of the cosmological magnetic field, and uq is the observed 
frequency of the radiation. Note that equation (|2|) is independent of cosmological parameters. 
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For a primordial field of 1CT 9 G which could result in the observed galactic field, we therefore 
expect a rotation measure of order 1.6 deg cm -2 = 280 rad m~ 2 . This rotation is considerable 
by astrophysical standards and could in principle be measured. The exact value of the 
rotation measure is, however, sensitive to the growth history of the microwave background 
polarization through the surface of last scatter. A larger Faraday rotation is expected in 
anisotropic cosmological models, previously investigated by Milaneschi & Fabbri (1985). 

In this paper, we perform a detailed calculation of the above Faraday rotation signal. 
We work in the context of inflation-type models with adiabatic initial fluctuations, and 
assume no early reionization. The rotation generated at the surface of last scatter depends 
only on the ionization history and is insensitive to most cosmological parameters. Section 
2 develops the formalism for calculating the microwave background polarization including 
Faraday rotation. In Section 3 we present numerical results, including the dependence of 
the rotation angle on the mean mass density and baryon density of the universe. Finally, 
Section 4 discusses the prospects for detecting this effect. 



The evolution of the cosmic microwave background is described by a set of radiative 
transfer equations for the Fourier modes of the radiation brightnesses A/(k, q, 77), Ag(k, q, 77), 
and A[/(k, q, 77), where the subscripts refer to the standard Stokes parameters, the Fourier 
mode wavevector is given by k, the radiation propagation direction is given by the unit vector 
q, and rj = J(l + z)dt is conformal time. For pure blackbody fluctuations, the temperature 
deviation is AT/T = Aj/4, with T = 2.726 ± 0.010 K (Mather et al. 1993) the mean 
temperature. In the case of no magnetic fields, the transport equations depend only on 
\i = cos(k • q), the cosine of the angle between a given Fourier mode and the propagation 
direction. Also, the isotropy of the Universe implies that the equations depend only on 
k — |k| and not on the direction of the Fourier component. In the following calculation, 
we assume a uniform magnetic field and extract mean results by averaging over the entire 
sky; then the equations still depend only on k and fi and not on q and k separately. This is 
equivalent to assuming that the magnetic field is coherent on the scale of the width of the 
last scattering surface, a comoving scale of about 5 Mpc. This assumption is natural if the 
coherent magnetic field observed in galaxies came from a primordial origin, since galaxies 
were assembled from a comoving scale of a few Mpc. 

Including the Faraday mixing term, the radiative transport equations in comoving 
coordinates are (Bond & Efstathiou 1984; Kosowsky 1996) 
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Formalism 



Aj + ifc//(Aj-4$) = -4^-fA 



./ - A /0 + 4v b fi - -P 2 (ll)(A I2 + A Q2 - Aqq) 
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A Q + ikfiAg 



A Q + - (1 - PM) (A I2 + Ag 2 — A Q0 ) 



2u; B A l/ 



Aiz + ifc/xA^ = -f A v - 2u B A Q , (3) 

where \1/ and $ are scalar metric perturbations in the Newtonian gauge, Vf, is the baryon 
velocity, and f is the differential optical depth, defined by f = x e n e a T a/ao with x e the 
ionization fraction, n e the total electron density, ot the Thomson cross-section, and a the 
scale factor normalized to ct today. Dots over quantities represent derivatives with respect to 
conformal time 77. The numerical subscripts on the radiation brightnesses indicate moments 
defined by an expansion of the directional dependence in Legendre polynomials Pe(f^): 

A It {k) = \ f 1 dnP t (n)&i(h,n). (4) 

2 J-l 

A detailed derivation of equations (f|) can be found in Kosowsky (1996). 

The mixing terms between Aq and Au account for the effect of Faraday rotation; these 
terms follow directly from the definition of the Stokes parameters. The conformal Faraday 
rotation rate is given by 

dtp dip a 

Wb = T = ' ( 5 ) 
dr] dt do 

where ip is the rotation angle of the polarization vector. For a given magnetic field and 
direction of photon propagation, the rotation rate dip/dt is given by equation ([!]). The 
time dependence of the gravitational potentials $ and ^ and the ionization fraction x e then 
completely determine the evolution of the radiation brightnesses. 

Hu and Sugiyama (1995a,b) have demonstrated how to solve the radiative transport 
equations semi-analytically, using analytic fits to the evolution of the potentials and the 
ionization fraction. This formalism has recently been extended to include polarization 
(Zaldarriaga & Harari 1995). Essentially, the tight-coupling solution describing the 
primordial plasma is modified by a damping factor to account for diffusion damping through 
recombination, and the resulting photon distribution free-streams to the present epoch. We 
use this approach to give the temperature fluctuations, which then source the polarization 
fluctuations. Here we present the modifications to the tight-coupling solution necessary 
to incorporate Faraday rotation; for other details of the calculation, see Hu & Sugiyama 
(1995a,b). 

The polarization brightness Aq is sourced by the function S p = A I2 + Aq 2 — Aq , and 
Ajj is generated as Aq and Ajj are rotated into each other. In the absence of magnetic fields, 
Ajj retains its tight-coupling value of zero. Expanding to second order in the tight-coupling 
parameter f _1 gives the evolution equation (Zaldarriaga & Harari 1995) 

3 2 

S p + —fSptt -ikAn, (6) 
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S p (k, V ) « like 3rM/1 ° [" drj' e - 3rW)/10 A I1 (k, V '). (7) 



with the solution 



5 Jo 
Note that the total optical depth r is defined as 

rv* 

t(v) = / Hv)dV, ( 8 ) 

with r]* the conformal time of recombination, giving dr/dr] = —f. We can rewrite the 
equations for the polarization brightnesses in a simple form with the change of variables 

A Q = e^~ T A Q , 
A v ee e^^Au, 

S p ee e ifcw - T ^ p , (9) 

giving the evolution equations 

3 

A Q (k,fx,rj) = --f(r])(l - fi 2 )S p (k,r]) + 2u B (fi,ri)A u (k, fi,rj), 

Au{k,fi,Ti) = -2aj B (n,r})A Q (k,n,ri). (10) 

We are interested in the case of small Faraday rotation; the first-order iterative solution to 
equations flip]) is 



Aq(* )A i,ij) = —(l-^jf'di/fCi/Je^^^C*;,!/), 

A,(M,^) = 3e3 ( B -q) (1 _ ») r dl / f (i/) fw^^i/l. (11) 

Equations (0) and flllD determine the polarization of the radiation as a function of frequency, 
given the ionization history x e and the temperature dipole brightness An. 

Using the polarization brightnesses in fc-space, we next obtain an expression for the 
polarization vector in x-space. Before performing the Fourier integral, the various spherical 
coordinate systems defined for each k-mode must be rotated to a common system since the 
definitions of the Stokes parameters Q and U depend on the orientation of the coordinate 
system (for details, see Kosowsky 1996). The real-space polarization fluctuations are given 
by 

Q(X ; M) = iEe ikX [A Q (k,^,0')cos2e / + A l/ (k,^,0')sin2e1, 
Jo 4 k 

U{ *' ' 0) = ^Ee lkX [-AQ(k,^,0')sin2e' + A [/ (k,^,0')cos2r], (12) 
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where (0', 0') represents the same direction as (0, 0) except in the coordinate system defined 
by the mode k. The angle £' is the rotation angle to align the coordinate systems; an 
expression for £' is given in Kosowsky (1996) but will not be used in what follows. Finally, 
the polarization vector P can be constructed from these quantities as 



P(M,. 



l 

72 



§y/P{P + Q)+4>lLy/P{P-Q) 



(13) 



with P = ^Q 2 + U 2 . The angle between two polarization vectors (pi 2 follows from 



2 (Pi • P 2 ) 

COS^ (pi2 = 



i + 



Q1Q2 + U1U2 
P1P2 



(14) 



The predictions of a cosmological model are only statistical in nature. For a given 
magnetic field, the observationally interesting quantity is the expectation value of the 
Faraday rotation angle. Calculationally, it is useful to consider only averages of quantities 
quadratic in the brightnesses given above; we thus consider the observable 



(cos cp 12 ) 



1 + 



(Q1Q2 + U1U2) 

(P 1 2 )l/2 (P 2 > l/2 



(15) 



The averages can be calculated explicitly in terms of the brightnesses by replacing the 
averages with the integral V^ 1 J dx. and replacing the Fourier sums with integrals, J2k - ► 
[V/(27r) 3 ] J dk, where V is a volume normalization factor. We also include an exponential 
beam suppression to account for a Gaussian beam of width a (Kolb & Turner 1990). The 
necessary averages can be written in terms of the two integrals 



/l /*oo 
dp J k 2 dkexp[-4:k 2 HQ 2 (l - A t 2 )cr 2 ]|A Q (A;,At)p 



/l re. 
dfi J o k 2 dkexp[-4k 2 Hv 2 (l - i2 2 )a 2 }\Au(k,i2)\ 2 



(16) 



with H the Hubble constant; a long calculation gives the final equations 

V 



(Q1Q2 + u ± u 2 ) 
( p l) 

(P2 2 ) 



64tt 2 
V 

64T 2 
V 

64T 2 



(Iq + flu) 
(Iq + Iu) 

(h + flu), 



(17) 



where the subscripts 1 and 2 refer to two different frequencies v\ < 1/2, all of the brightnesses 
are evaluated at vi, and / = v\jv\. Since we are comparing the polarization at two 
frequencies in the same direction on the sky, the formulas look like the temperature 
correlation function at zero separation. 
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In the perturbative limit of small rotation angles, equation (15) reduces to 

<^ 2 > 1/2 - ~ 

Equations ([TT|), (|T6|), and (0) imply that the rms rotation angle is proportional to the 
magnetic field and inversely proportional to the square of the frequency, as expected. An 
average over all directions of observation will reduce the above estimate by a factor of a/2, 
due to the changing orientation of the magnetic field. 



3. Results 



We solve the equations in the previous section numerically to determine l\j / Iq and thus 
the proportionality factor between the mean Faraday rotation signal and the quantity Bjv 2 . 
Normally, the radiative transfer equations for the microwave background are expanded in 
a moment hierarchy to eliminate the fi dependence; however, in this case, the additional 
directional dependence of the magnetic field direction complicates the moment expansion. 



We perform the integrals in equations (pTTj) on a (/i,fc) grid. The necessary cosmological 
inputs are the dipole component of the radiation intensity, A/i, and the differential optical 
depth through recombination, f. 

We choose as an illustrative model "standard" cold dark matter, but the nature of 
the perturbations has little effect on our results; the size of the radiation dipole and the 
details of recombination affect the final answer, but neither of these depends strongly on the 
cosmological model. The matter density Qoh 2 and the baryon density Qbh 2 have a mild effect 
on recombination, which we include through numerical evolution of the free electron density. 
Hu and Sugiyama (1995a) give an analytic approximation for An which we use here. We 
evolve the electron density using a numerical code incorporating the recombination physics 
detailed in Hu et al. (1995). 

Figure 1 displays the evolution of the polarization brightnesses Aq and Au through 
the last scattering surface for \x = 0.5 along with the differential visibility function fe~ T . 
At early times, the tight coupling between the photons and the baryons prevents the 
development of polarization. As decoupling proceeds, the photons begin to free-stream, 
generating a quadrupole perturbation which sources the polarization. However, the induced 
Faraday rotation depends on the free electron density, which drops to negligible values as 
recombination ends. Rotation is generated during the brief period of time when the free 
electron density has dropped enough to end tight coupling but not so much that Faraday 
rotation ceases. Figure 1 shows that the generation of rotation lags behind the generation 
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Fig. 1. — The evolution of the polarization brightnesses, for k = 0.16 Mpc~ and \i = 0.5 (in 
arbitrary units). Also plotted as a dotted line is the differential visibility function fe~ T in 
units of Mpc -1 . 

of polarization, and that rotation is essentially completed by z = 1100 whereas polarization 
generation continues almost until z = 900. 

Figure 2 displays Iu/Iq m equation fll8[) as a function of the cosmological parameters 
Q h 2 , the mass density parameter, and Q^h 2 , the baryon density parameter, for B = 10~ 9 
G, V\ = 30 GHz, and a beam-width of 0.5°. The beam-width only affects the signal-to-noise 
ratio of the polarization signal and has virtually no effect on the size of the rotation angle. 
As expected in the Introduction, the value of (ip 2 ) oc Iu/Iq is nearly independent of these 
cosmological parameters. To within a 10% correction for the effects of Qbh 2 and VLqIi 2 (in 
the range VL^h 2 > 0.007 and Q h 2 < 0.3), the rotation angle is 

in good agreement with equation (^). We have included here a factor of l/y/2 due to an 
average over all orientations of the magnetic field. 
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Fig. 2. — The ratio Iq/Iu as a function of £l h 2 and fi b /i 2 , for a magnetic field _B = 10 9 
G and an observed frequency of 30 GHz. The dotted lines from bottom to top are for 
Q b h 2 = 0.005, 0.0075, and 0.01. The solid lines from top to bottom are for Q b h 2 = 0.0125, 
0.015, 0.02, and 0.025. 



We have calculated the rotation of the microwave background polarization vector due to 
a primordial magnetic field at the surface of last scatter. The magnetic field was assumed to 
be coherent on the width of the last scattering surface, corresponding to a comoving scale of a 
few Mpc, which is conveniently the length scale from where galaxies are assembled. We have 
found that the mean Faraday rotation has a size of about 1° for a cosmological magnetic 
field of Bo = 10~ 9 G and an observed frequency of 30 GHz. An optimized experiment 
could measure the orientation of the polarization vector at one frequency which is as low as 
practicable, and at a second, somewhat higher frequency. If the frequency ratio is 2, then 
75% of the total rotation is observed between the two frequencies. While extracting this 
rotation in a single observation direction is unlikely, a statistical detection averaged over 
many observation directions is possible. 



4. 



Discussion 
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The peak amplitude of the expected polarization fluctuations is an order of magnitude 
lower than that of the microwave temperature fluctuations, of order 10~ 6 (see, e.g., 
Crittenden, Davis, & Steinhardt 1993; Frewin, Polnarev, & Coles 1994). To measure a 
Faraday rotation of 1° requires another factor of 10 2 in sensitivity. The total sensitivity to 
the Faraday rotation signal is proportional to the raw pixel sensitivity and to the square 
root of the number of pixels. Currently envisioned satellite experiments are designed with 
the primary goal of mapping the microwave sky temperature down to small angular scales; 
given amplifiers of a certain sensitivity, it is advantageous to design an experiment with a 
signal-to-noise ratio per pixel of around unity and as many pixels on the sky as possible 
(Knox 1995). Such a mapping experiment would require of order 10 6 pixels to detect a 
rotation angle of 1°. This corresponds to a beam size of order 10' which is at the limit 
of current design proposals. It is not unreasonable to expect that in the coming years, raw 
sensitivity will improve so that systematic effects become the dominant obstacle to detecting 
the Faraday rotation signal. As with temperature measurements, the ultimate limitation will 
be contamination from foreground sources. While little is known about polarization sources 
at microwave frequencies, estimates suggest that in the frequency window of 30-80 GHz 
the signal should be dominated by the microwave background contribution (Timbie 1995). 
Experiments which focus on high signal-to-noise observations of small patches of the sky may 
also prove useful; their feasibility could be first demonstrated by searching for a foreground 
rotation signal in the direction of galaxy clusters with detected rotation measures from radio 
source observations (see, e.g. Dreher et al. 1987; Perley & Taylor 1991; Taylor & Perley 1993; 
Ge & Owen 1993). In fact, microwave background polarization measurements at multiple 
frequencies can potentially map the magnetic field distribution in cluster environments where 
the rotation measure is already known to be large (Loeb & Kosowsky 1996). 

Will foreground contamination of the rotation measure be a considerable obstacle to 
measuring the Faraday rotation signal from primordial recombination? If the universe 
was reionized at a redshift z <C 100 x f2 1/,3 (f2b/i/0.02)~ 2 / 3 , then the optical depth through 
reionization is much smaller than unity and the additional Faraday rotation contribution is 
small. The contribution from late reionization would be further reduced if the coherence 
length of the primordial field is much smaller than a Gpc. In fact, the cumulative rotation 
measure of the intergalactic medium out to z ~ 3 is known to be smaller than a few 
rad m~ 2 from samples of high-redshift objects (see Valee 1990, Kronberg 1994, and references 
therein). In addition, the rotation measure of the Milky- Way galaxy was measured at high 
latitudes and found to be ^ 20 rad m~ 2 (Spitzer 1978; Simard-Normandin & Kronberg 1980). 



According to equation (|19l) , these foregrounds should allow the detection of a cosmological 
field as small as 10~ 10 G. 
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